4,876 research outputs found
A high-order approach to elliptic multiscale problems with general unstructured coefficients
We propose a multiscale approach for an elliptic multiscale setting with
general unstructured diffusion coefficients that is able to achieve high-order
convergence rates with respect to the mesh parameter and the polynomial degree.
The method allows for suitable localization and does not rely on additional
regularity assumptions on the domain, the diffusion coefficient, or the exact
(weak) solution as typically required for high-order approaches. Rigorous a
priori error estimates are presented with respect to the involved
discretization parameters, and the interplay between these parameters as well
as the performance of the method are studied numerically
Multiscale scattering in nonlinear Kerr-type media
We propose a multiscale approach for a nonlinear Helmholtz problem with possible oscillations in the Kerr coefficient, the refractive index, and the diffusion coefficient. The method does not rely on structural assumptions on the coefficients and combines the multiscale technique known as Localized Orthogonal Decomposition with an adaptive iterative approximation of the nonlinearity. We rigorously analyze the method in terms of well-posedness and convergence properties based on suitable assumptions on the initial data and the discretization parameters. Numerical examples illustrate the theoretical error estimates and underline the practicability of the approach
Computational multiscale methods for linear heterogeneous poroelasticity
We consider a strongly heterogeneous medium saturated by an incompressible
viscous fluid as it appears in geomechanical modeling. This poroelasticity
problem suffers from rapidly oscillating material parameters, which calls for a
thorough numerical treatment. In this paper, we propose a method based on the
local orthogonal decomposition technique and motivated by a similar approach
used for linear thermoelasticity. Therein, local corrector problems are
constructed in line with the static equations, whereas we propose to consider
the full system. This allows to benefit from the given saddle point structure
and results in two decoupled corrector problems for the displacement and the
pressure. We prove the optimal first-order convergence of this method and
verify the result by numerical experiments
Localized implicit time stepping for the wave equation
This work proposes a discretization of the acoustic wave equation with
possibly oscillatory coefficients based on a superposition of discrete
solutions to spatially localized subproblems computed with an implicit time
discretization. Based on exponentially decaying entries of the global system
matrices and an appropriate partition of unity, it is proved that the
superposition of localized solutions is appropriately close to the solution of
the (global) implicit scheme. It is thereby justified that the localized (and
especially parallel) computation on multiple overlapping subdomains is
reasonable. Moreover, a re-start is introduced after a certain amount of time
steps to maintain a moderate overlap of the subdomains. Overall, the approach
may be understood as a domain decomposition strategy (in space and time) that
completely avoids inner iterations. Numerical examples are presented
Computational Multiscale Methods for Linear Poroelasticity with High Contrast
In this work, we employ the Constraint Energy Minimizing Generalized
Multiscale Finite Element Method (CEM-GMsFEM) to solve the problem of linear
heterogeneous poroelasticity with coefficients of high contrast. The proposed
method makes use of the idea of energy minimization with suitable constraints
in order to generate efficient basis functions for the displacement and the
pressure. These basis functions are constructed by solving a class of local
auxiliary optimization problems based on eigenfunctions containing local
information on the heterogeneity. Techniques of oversampling are adapted to
enhance the computational performance. Convergence of first order is shown and
illustrated by a number of numerical tests.Comment: 14 pages, 9 figure
Reconstruction of quasi-local numerical effective models from low-resolution measurements
We consider the inverse problem of reconstructing an effective model for a prototypical diffusion process in strongly heterogeneous media based on low-resolution measurements. We rely on recent quasi-local numerical effective models that, in contrast to conventional homogenized models, are provably reliable beyond periodicity assumptions and scale separation. The goal of this work is to show that the identification of the matrix representation of these effective models is possible. Algorithmic aspects of the inversion procedure and its performance are illustrated in a series of numerical experiments
An improved high-order method for elliptic multiscale problems
In this work, we propose a high-order multiscale method for an elliptic model
problem with rough and possibly highly oscillatory coefficients. Convergence
rates of higher order are obtained using the regularity of the right-hand side
only. Hence, no restrictive assumptions on the coefficient, the domain, or the
exact solution are required. In the spirit of the Localized Orthogonal
Decomposition, the method constructs coarse problem-adapted ansatz spaces by
solving auxiliary problems on local subdomains. More precisely, our approach is
based on the strategy presented by Maier [SIAM J. Numer. Anal. 59(2), 2021].
The unique selling point of the proposed method is an improved localization
strategy curing the effect of deteriorating errors with respect to the mesh
size when the local subdomains are not large enough. We present a rigorous a
priori error analysis and demonstrate the performance of the method in a series
of numerical experiments.Comment: 18 pages, 4 figure
Reconstruction of Quasi-Local Numerical Effective Models from Low-Resolution Measurements
We consider the inverse problem of reconstructing an effective model for a prototypical diffusion process in strongly heterogeneous media based on coarse measurements. The approach is motivated by quasi-local numerical effective forward models that are provably reliable beyond periodicity assumptions and scale separation. The goal of this work is to show that an identification of the matrix representation related to these effective models is possible. On the one hand, this provides a reasonable surrogate in cases where a direct reconstruction is unfeasible due to a mismatch between the coarse data scale and the microscopic quantities to be reconstructed. On the other hand, the approach allows us to investigate the requirement for a certain non-locality in the context of numerical homogenization. Algorithmic aspects of the inversion procedure and its performance are illustrated in a series of numerical experiments
- …